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ABSTRACT 


The  Linearly  Tapered  Slot  Antenna  (LTSA)  had  been 
investigated  and  developed  experimentally;  its  applications 
have  primarily  been  based  on  empirical  designs.  An  accurate 
theoretical  model  based  on  Moment  Methods  (MM)  is  developed 
here  to  study  the  radiation  characteristics  of  the  LTSA. 
Using  the  MM  solutions  to  the  reaction  integral  equation, 
this  thesis  presents  an  analysis  to  model  and  explain  the 
LTSA  behavior.  The  effects  of  variable  design  parameters  on 
radiation  patterns  are  studied.  Discussion  is  augmented  by 
relating  predicted  radiation  patterns  to  calculated  current 
distributions  on  the  antenna  surface.  Conclusions  are  made 
regarding  optimum  designs  for  the  LTSA.  Relevant 
observations  are  made  concerning  the  extensive  computational 
tasks  and  the  computer  resources  required  for  the  MM  model. 
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I. 


INTRODUCTION 


The  many  practical  advantages  of  the  Linearly  Tapered 
Slot  Antenna  (LTSA)  are  based  on  its  low  cost,  light  weight 
and  relative  simplicity.  Specific  designs  of  LTSA  arrays 
have  been  implemented  for  millimeter  wave  imaging  systems  and 
more  applications  are  inevitable.  However,  since  the 
refinement  of  an  accurate  theoretical  model  has  lagged  its 
development,  LTSA  uses  have  been  based  on  empirical  designs. 

Based  on  the  proven  Moment  Methods  (MM)  model,  developed 
by  Janaswamy  [Ref.  1],  this  paper  presents  an  analysis  of  the 
LTSA  radiation  pattern.  The  model  permits  a  parametric  study 
of  the  antenna.  The  physical  antenna  variables  used  in  the 
LTSA  modeling  are  illustrated  in  Figure  1.1. 

Since  the  method  of  this  research  is  the  MM  model,  a 
general  review  of  MM  and  a  specific  discussion  of  the  LTSA 
model  are  included  in  Chapters  2  and  3  respectively.  It  is 
intended  that  these  chapters  address  both  the  accuracy  of  MM 
as  well  as  the  limitations,  particularly  with  respect  to 
solution  by  computer.  By  way  of  examples,  Chapter  2  is 
demonstrative  of  the  solution  by  the  Moment  Method.  The  MM 
concept  is  then  advanced  in  Chapter  3  to  explain  how  the 
reaction  integral  equation  is  solved  in  the  LTSA  model. 
Chapter  3  further  serves  to  discuss  the  model  in  terms  of 
computer  tasking  considerations  and  approximations  which 
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alleviate  computational  burden  with  little  or  no  sacrifice  in 
accuracy. 

Chapter  4  focuses  on  the  LTSA  modeling  results  and 
analyzes  the  observed  effects  of  changing  beamwidths  and 
sidelobe  suppression  in  the  E-  and  H-planes.  Of  particular 
interest  are  the  pattern  changes  which  result  from  various 
combinations  of  antenna  length  and  height,  but  effects  of 
changing  other  design  parameters  are  likewise  investigated. 
Corresponding  to  the  radiation  pattern  changes,  surface 
current  distributions  are  used  to  explain  some  phenomena 
which  appear  to  synergistically  create  the  different 
patterns.  (Appendix  C  contains  current  distribution  graphs 
for  the  respective  radiation  patterns  illustrated  in  Appendix 
A)  . 

Conclusions  are  summarized  in  Chapter  5.  The  thrust 
therein  is  to  provide  a  credible  spectrum  of  design 
considerations  for  the  LTSA.  It  is  fully  expected  that  this 
work  has  provided  the  necessary  data  and  rationale  for 
relevant  interpolation  as  well  as  extrapolation.  Moreover, 
the  methods  used  in  this  study  may  well  serve  similar  studies 
of  other  antennas  and  scatterers. 
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II.  THE  MOMENT  METHOD 


A.  CONCEPT 

Although  the  Moment  Method  (MM)  is  a  widely  used 
technique  in  solving  electromagnetics  problems,  only  the 
advances  in  computer  technology  of  recent  years  enable  the 
cost-effective  use  of  MM  to  pursue  such  rigorous  solutions. 
Specifically,  a  linear  operator  equation  with  an  unknown 
function  is  solved  by  expanding  the  unknown  function  and 
casting  it  into  a  system  of  linear  equations  which  is  then 
solved  by  known  techniques.  The  thrust  of  MM  is  to  maintain 
the  exactness  of  the  original  equation  and  to  solve  it 
numerically  or  by  other  means.  This  is  normally  done  at  the 
expense  of  additional  computer  central  processing  unit  (CPU) 
time.  However,  MM  generally  provides  a  more  accurate 
solution  in  comparison  to  other  techniques  which  reduce  and 
approximate  equations  solely  to  minimize  the  computational 
effort . 

The  prescription  for  solving  a  problem  by  MM  is  best 
illustrated  by  an  example.  Two  such  examples  are  chosen  for 
this  purpose.  The  first  problem  is  to  solve  for  f(x)  in  the 
Fredhom  integral  equation  of  the  first  kind 

1 

j  f ( x ) cos ( xy ) dx  =  (1-cos  y)/y2.  (2.1) 

J0 


4 


The  exact  solution  to  the  above  equation  is 


f(x)  =  1-x,  0<x< 1 .  (2.2) 

While  exploring  this  elementary  problem  in  detail,  it  is 
intended  to  demonstrate  the  MM  technique  and  to  highlight 
both  its  potential  advantages  as  well  as  its  shortcomings. 
The  second  sample  problem  is  to  solve  for  the  current 
distribution,  I(z')  on  a  short  dipole  using  Pocklington ' s 
integral  equation  [Ref.  2:p.  307]: 


1 

]  u)£o 


PL/2 


-L/2 


+  £2 


-E*(z) . (2.3) 


The  development  and  analysis  of  this  equation  can  be  found  in 
most  antenna  text  books  and  only  the  application  of  MM  is  of 
current  interest. 


B.  SOLUTION  OF  THE  FREDHOLM  INTEGRAL  EQUATION 

In  solving  the  Fredholm  integral  equation,  the  analytical 
approach  is  somewhat  awkward.  However,  MM  is  used  with 
principles  of  linearity  and  superposition  to  reduce  the 
functional  equation  to  a  matrix  equation  of  the  form  [A][x]  = 
[B]  [Ref.  3].  Then  known  techniques  are  used  to  determine  a 
unique  solution  for  [x]  =  [A]-1[B].  Equation  2.1  in  terms  of 
a  linear  operator,  L,  is  written  as 


where , 


1 

L  =  f  cos(xy)dx  (2.5) 

J0 

g  =  (1-cos  y)/y2.  (2.6) 

The  first  step  is  to  expand  f  into  a  set  of  basis 
functions,  <fi,  f2...fM}.  N  determines  the  size  of  matrix  to 
be  used  and  implies  both  the  degree  of  accuracy  and  the 
required  computational  time.  The  selection  of  fn  is  more 
crucial  and  should  be  selected  to  approximate  the  expected 
behavior  of  f(x).  Figure  2.1  shows  typical  functions  with 
which  to  represent  f,  where  N  =  5  and 

N 

f  =  Z  anf„ .  (2.7) 

n=l 

The  problem  reduces  to  solving  for  "best  fit"  coefficients, 
an  for  each  fn.  Substituting  equation  (2.7)  into  (2.4) 
yields 

N 

Z  a„«L(f„)  =  g.  (2.8) 

n=l 

There  are  now  essentially  N  unknowns,  and  the  number  of 
linearly  independent  equations  must  also  equal  N  for  a  unique 
solution.  This  is  achieved  by  selecting  a  set  of  weighting 
functions,  (wt,  wa...wM}  in  the  range  of  L  and  equating  the 
weighted  averages  of  equation  (2.8)  for  each  w„ , 
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1  1 

0  xl 

1  | 

x2  x3 

(a)  Dirac  delta  function 

1 

x4 

- i 

x5 

h 

P(x-xl) 

] _ _ 

1 

1  1 

1 

1 

0  xl 

x2  x3 

x4 

x5 

(b)  Pulse  function 


1  I  I  I  I  I 

0  xl  x2  x3  x4  x5 


(c)  Step  approximation 


0  xl  x2  x3  x4  x5 


(d)  Function  represented  by  triangles 


0  xl  x2  x3  x4  x5 

(e)  Piecewise  sinusoid  function 


Figure  2.1  Typical  Basis  Functions. 
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N 

£  an*<w„,  L  (  f  n )  >  =  <wm ,  g> .  (2.9) 

n=l 

Here  <  >  denotes  the  inner  product,  defined  as 

1 

<f,g>  =  f  f(x)  g(x)  dx.  (2.10) 

J  0 

For  M  equations  and  N  unknowns,  equations  (2.9)  is 
expressed  as 


[  1  mn  ]  [  OCn  ]  = 

[gmJ 

(2, 

.11) 

<Wi , 

Lf  i 

L  > 

<Wj, 

Lf 2> . . . <Wi , 

LfN> 

<Wj, 

Lf , 

L> 

<Ws, 

Lf  2^  .  .  • 

(2. 

12) 

<WM, 

Lf, 

L'  i  •  « 

— 

tr* 

t-h 

Z 

V 

1 _ 

[a 

o]  = 

=  [a. 

l  ,  CLz 

.  .  .  CLn  ]T 

(2. 

.13) 

t  Qrn 

]  = 

[9i  i 

g2 . 

•  • 5m  j T . 

(2. 

14) 

Thus,  the  unknown  coefficients  are  determined  by 

[an]  =  [lmn]-xtg«n].  (2.15) 

A  point  not  to  casually  dismiss  is  the  selection  of 
weighting  functions.  It  is  essential  that  both  fn  and  wn  be 
linearly  independent  sets  and  that  the  resulting  "1"  matrix 
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have  an  inverse  and  be  well  conditioned.  The  infinite 
possibilities  for  fr,  and  wn  allow  considerable  latitude  to 
accommodate  (a)  the  degree  of  accuracy  desired,  (b)  ease  of 
computation  and  (c)  an  acceptable  size  matrix  for  inversion. 
Choices  for  wn  should  be  generally  based  on  these  criteria 
while  selections  of  f„  are  to  help  approximate  the  expected 
or  assumed  behavior  f(x). 

For  illustrative  purposes,  the  above  Fredholm  integral 
equation  is  solved  with  N  pulse  functions  as  basis  functions. 
The  pulse  function  P(x)  is  defined  as 

(1  | x |  <  1/2N 

P(x)  =  J  .  (2.16) 

(0  | x |  >  1/2N 

Thus,  a  linear  combination  of 

f„  =  P(X-Xn)  (2.17) 

gives  a  step  approximation  to  f(x). 

Similarly,  N  weighting  functions  are  chosen  as  Dirac 
delta  functions: 

wm  =  6(x-Xm) ,  (2.18) 

with  xm  taken  as  the  midpoint  of  each  subinterval.  This 
choice  of  weighting  functions  is  also  known  as  the  point¬ 
matching  technique.  Physically  this  choice  implies  that  the 
original  integral  equation  (2.1)  is  now  imposed  at  the 
discrete  points  xm  (the  locations  of  the  delta  functions), 
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These 


rather  than  in  the  continuous  interval,  0<x<l. 
selections  serve  primarily  to  easily  demonstrate  the 

computations  rather  than  to  faithfully  reproduce  f(x). 

Figure  2.2  shows  a  pictorial  representation  of  the  basis  and 
weighting  functions. 

For  a  manageable  demonstration,  set  N=M=3.  The  Fredholm 
integral  equation  (2.1)  is  expressed  in  the  form  of  equation 
(2.9): 


3 

£  a„ • <6 ( y-xm) 
n=l 


%  a 

P(x-Xn)  cos(xy)dx> 

O 


=  <6(y-xm),  ( l-cos/ya> 


(2.19) 


where  the  inner  product  is  computed  using  dummy  variable,  y. 
Therefore,  from  equation  (2.12),  each  matrix  element 


1  1 

lmn  =  I  6(y-xm)  |  P(x-x„)  cox(xy)dxdy 

J  o  J  0 


n/3 

(n-l)/3 


cos (xy )dxdy 


n/3 

(n-l)/3 


cos(x*xm)dx 


f 
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1 

x„ 


nxm  (n-l)xr 

(sin  -  -  sin  - 

3  3 


(2.20) 


In  the  same  manner, 


Qm  = 


6(y-xm) 


1  -  cos  y 


1  -  cos  x„ 


dy  = 


(2.21) 


Solving  equations  (2.20),  (2.21)  and  (2.15)  yields 


.3332 

.3321 

.3301 

[lmn] 

= 

.3318 

.3226 

.3045 

(2.22) 

.3291 

.3038 

.2553 

[gm]  = 

[ . 

.4988 

.4897 

.4717  ]T 

(2.23) 

[a„]  = 

[  .8319 

.5647 

.  1034  ]T. 

(2.24) 

This  solution  is  depicted  with  the  exact  solution  (f(x)  = 
1-x)  in  Figure  2.3. 

Figure  2.3  seems  to  suggest  that  the  accuracy  of  the 
numerical  solution  may  be  improved  by  increasing  the  number 
of  basis  functions.  While  this  may  be  true  in  general,  it 
needs  a  high  degree  of  precision  to  actually  achieve  it.  For 
a  given  N,  the  percent  error  is  less  if  the  precision  is 
increased.  Conversely,  for  a  given  precision,  increasing  N 
eventually  results  in  an  increase  in  percent  error.  The 
accuracy  will  also  be  affected  by  moving  the  matching  point, 
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-  Numerical  solution 


z 


Exact  solution 


l 


x 


)  and  Exact  Solutions  for  the 


i.e.,  by  shifting  the  weighting  functions  6(x-xm)  within  the 
corresponding  subinterval 


m-1  m 

-  <  xm  <  (2.25) 

N  N 


In  another  approach  known  as  Galerkin's  technique,  the 
set  of  weighting  functions  is  the  same  as  the  set  of  basis 
functions,  i.e.,  w„  =  f„.  This  results  in  a  symmetric  NxN 
matrix  and  reduces  the  required  number  of  computations  needed 
for  filling  the  matrix. 

Reviewing  a  few  different  MM  solutions  to  the  Fredholm 
integral  provides  greater  insight.  Various  MM  approaches  and 
results  are  itemized  in  Table  2.1.  For  this  problem  and  the 
limited  trials  indicated,  a  few  conclusions  are  made: 

*  The  noticeably  low  value  of  N  is  the  largest  which  can 
be  used  before  the  matrix  becomes  ill-conditioned.  This 
is  attributable  to  the  available  limits  of  precision  in 
the  mantissa. 

*  Use  of  Galerkin's  technique  does  not  appear  to  have  a 
dominant  effect  on  accuracy. 

*  Selection  of  apparently  suitable  basis  functions  does 
not  guarantee  convergence,  since  any  N  too  large  may 
result  in  an  ill-conditioned  matrix. 

*  Shifting  the  weighting  functions  in  the  point  matching 
technique  does  not  significantly  affect  the  solution  for 
la„]. 

The  ambiguous  tone  of  these  conclusions  alludes  to  the 
obvious  point  that  few  solutions  offer  the  neatness  of  f(x)  = 
1 — x .  What  is  to  be  gleaned  from  this  example  is  that  there 
are  limitations  to  the  MM  process.  No  exact  guidelines  exist 
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TABLE  2.1  REPRESENTATIVE  SELECTIONS  OF  BASIS  AND 

WEIGHTING  FUNCTIONS  TO  SOLVE  THE  FREDHOLM 
INTEGRAL. 


Basis 

Function 

Weighting 

Function 

N 

Precision 

Solution[aI,]‘r 

Pulse 

Dirac  Delta  enforced 

3 

Single 

[.8319  .5647 

§  center  of  sub- 

.  1034  ] 

interval 

Pulse 

Dirac  Delta  enforced 

3 

Single 

[.8343  .5611 

@  left  of  subinterval 

. 1046  J 

Pulse 

Dirac  Delta  enforced 

5 

Double 

[.9622  .6283 

@  left  of  subinterval 

.5167  .3288 

.0641] 

Pulse 

Pulse  Galerkin 

3 

Single 

[.8449  .5451 

. 1101  ] 

Pulse 

Pulse  Galerkin 

3 

Double 

[.8379  .5557 

.1064] 

Powers  of  x 

Dirac  Delta  enforced 

3 

Single 

[  .9778 

(x^x1 .  .  . 

@  xm= ( 2m-l ) /N 

-.8841 

xM-') 

.1072] 

Powers  of  x 

Dirac  Delta  enforced 

2 

Single 

[  .9991 

( x° , X1  .  .  . 

@  xm= ( 2m-l ) /N 

-.9881] 

X”-1) 

Powers  of  x 

Powers  of  x(Galerkin) 

2 

Single 

[  .9974 

(x° ,  X1 . . . 

@  xm= ( 2m- 1 ) /N 

-.9946] 

xM~l) 
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for  the  best  choice  of  expansion  and  weighting  functions  for 
a  given  integral  equation.  In  spite  of  this  shortcoming,  MM 
is  much  more  useful  in  solving  an  integral  equation, 
particularly  when  it  cannot  be  solved  in  a  closed-form,  e.g., 
A^cos  (wt).  This  point  is  reinforced  in  the  next  example. 

C.  SOLVING  FOR  CURRENT  DISTRIBUTION  ON  A  DIPOLE 

Consider  now  an  electromagnetics  application  of  MM  to 
solve  for  current  distribution  on  a  linear  dipole.  The 
motivation  for  finding  current  distribution  on  any  antenna  is 
to  eventually  solve  for  the  radiation  pattern  and/or  the 
input  impedance.  Since  this  is  the  idea  of  the  LTSA  model 
specifically  discussed  in  Chapter  3,  an  example  of  MM  to 
solve  the  Pocklington  integral  equation  is  worthwhile. 

Using  the  conventional  assumptions  for  thin  wire  antennas 
and  illustrations  of  Figure  2.4,  the  problem  statement  is  to 
solve  for  the  current,  I(z')  on  a  short  dipole  with  length,  L 
=  .lA  and  radius,  a  =  .005  A,.  Assume  a  1  volt  center-fed 
excitation . 

In  a  simpler  form,  equation  (2.2)  is  written  as 

L/2  , 

I(Z'  )  K  (  z  ,  z  '  )  dz’  =  -E't  (2.26) 

L/2  Z 

where  the  kernel, 

K(  Z  ,  Z  '  )  =  -  +  S2^(  z  ,  Z  '  )  (2.27) 

<?Z2 


-0_  f 

]w£o  J 
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z 


Figure  2.4  Linear  Dipole  Illustrated  as  a  Wire  With  Surface 
Current  Density  J  s  and  Observation  Point  on  the  Surface. 
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and  the  free  space  Green's  function  is 

e  -  3  ram 

^(Z,Z‘  )  =  - .  (2.28) 

4rtR 

Equation  (2.25)  is  expanded  and  simplified  for  integration: 

0-  3 

K( z , z  1  )  =  - [(1  +  j  GR ) *  ( 2R2  -  3a2)  +  azfi2R2 ] . (2.29) 

4rRs 

It  is  readily  apparent  that  the  required  complex  integration 
is  best  done  numerically  and  is  well-suited  for  a  modest 
computer  program.  However,  for  present  purposes  only,  the  MM 
concept  and  results  need  be  explained. 

Following  the  same  procedure  as  in  the  previous  example, 
the  steps  are  simple:  (a)  expand  I(z')  in  a  set  of  basis 

functions;  (b)  obtain  linear  equations  by  weighting  functions 
and  computing  inner  products;  and  (c)  use  matrix  algebra  to 
solve  for  the  unknown  coefficients  of  the  basis  functions. 

For  simplicity  sake,  N=5  is  chosen  to  limit  the  matrix 
dimensions.  For  interpretive  convenience,  I(z')  is 

approximated  by  a  series  of  pulse  functions  (see  Figure  2.5): 

5 

I(Z’)  =  S  InFn(Z')  (2.30) 

n=l 

where  i„  is  a  complex  unknown  number  and 

II  for  z'  in  z'n 

(2.31) 

0  elsewhere. 
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proximation  to  an  Actual  Current 


9 


Substituting  equation  (2.30)  into  (2.26)  yields 

1  5  \^/7  i 

-  2  I„  Fn(z' )K(Z,Z’ )dz’  ~  -  E  (2).  (2.32) 

jw£o  n=l  ,  —  L./2  z 

Selection  of  weighting  functions  is  influenced  by  the  desire 
to  make  computations  simple.  Therefore,  the  point-matching 
technique  is  used.  The  weighting  functions  are 

wm(  z )  =  5 ( z-zm) .  (2.33) 

Taking  the  inner  products  in  equation  (2.32)  produces 

1  5  fL'2  i 

-  2  I„  F„(Z'  )K(zm,Z'  )dz'  *  -  E  ( Z m )  •  (2.34) 

jw£o  n=l  J-1,/2  z 

The  subscript,  m,  is  considered  the  match  point  index  because 
it  is  associated  with  the  observation  point  at  which  the  m*'” 
equation  is  valid.  The  source  point  index,  n,  is  so  called 
because  it  refers  to  the  field  from  the  current  on  the  nr’’ 
segment.  For  mathematical  convenience,  the  locations  of 
enforcement  are  at  the  centers  of  each  segment.  (It  is 
important  to  recognize  that  while  mathematically  simple, 
point  matching  actually  relaxes  the  boundary  conditions 

i  s 

E  +  E  =  0  (2.35) 

z  z 
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between  selected  points.  Such  simplification  should  be 
questioned  before  hastily  using  this  technique  in  other 
antenna  problems). 

Equation  (2.34)  is  solved  as  a  matrix  equation: 

[Zmn][I„]  =  [Vm]  (2.36) 


where 


Ju)£o  J 


nL/  2 


(  n-  1  >  Z>/2 


Fn(z 1 ) K( zm, z ' ) dz ' 


(2.37) 


i 

Vm  =  -E  ( z  m ) •  (2.38) 

z 


Note  the  mathematical  convenience  of  (2.36)  which  is 
analogous  (but  not  equivalent)  to  Kirchhoff's  network 
equations.  Because  of  the  analogy,  however,  [Zm„], [In]  and 
ivm]  are  respectively  referred  to  as  generalized  impedance, 
current  and  voltage  matrices. 

Before  solving  equation  (2.36),  [VmJ  must  be  defined  for 
the  given  problem.  Source  modeling  is  a  subject  worthy  of 
exploration  beyond  this  study  [Ref.  2].  Two  most  popular 
versions  are  (a)  the  delta  gap  generator  that  corresponds  to 
feeding  by  an  ideal  voltage  source  and  (b)  a  magnetic  frill 
generator  that  corresponds  to  feeding  by  a  coax  cable  from  a 
ground  plane.  Here,  a  1  volt  center  fed  excitation,  when 
modeled  with  a  frill  generator  and  weighted  with  delta 
testing  functions,  results  in 
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I  V„  ]  = 


0 . 502 
3 . 248 
70.55 
3 . 248 
.  502 


z.  -179.69° 
2:  -179.95° 
-179.99° 
z.  -179.95° 
z.  -179.69° 


(2.39) 


The  above  choice  of  uniform  segmentation  of  the  wire  results 
in  an  impedance  matrix  that  is  Toeplitz  symmetric.  (In  a 
Toeplitz  matrix,  Zmn  =  Zi, im-rn+t  for  m£2,  nZl.)  Thus,  the 
first  column  (identical  to  the  first  row)  is  all  that  is 
needed  to  fill  Zmn. 


679 . 2 

2i  90.01° 

292.4 

2i  -90.03° 

33 . 01 

z.  -90.27° 

9.74 

/.  -90.01° 

4.24 

21 -92.08° 

Solving  (2.36)  for  [In]  yields 


.81 

21  89.54° 

1.54 

2i  89.64° 

[I~] 

=  103  • 

2.44 

z.  89 .75° 

1 . 54 

z.  89 . 64° 

.81 

z.  89,54° 

solution 

for  [InJ 

is 

consistent 

(2.41) 


with  expected 

results.  The  current  distribution  is  generally  triangular. 
The  input  impedance  is  largely  capacitive  with  a  real  part 
less  than  2  ohms.  Nevertheless,  a  more  accurate  solution  for 
[In]  is  obtainable  with  a  larger  N.  Figure  2.6  shows  how  a 
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L IX 


Figure  2.6  Current  Distributions  on  a  Half-wave  Dipole  for  Various 
Numbers  of  Basis  Functions  and  a  Unit  Volt  Source  [Ref.  2:  p.314]. 
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solution  for  [ln]  converges  by  increasing  N  in  a  similar 
problem  [Ref.  2:p.  314]. 

The  purpose  of  this  example  has  been  to  demonstrate  the 
MM  application  to  a  simple  antenna  problem.  It  is  necessary 
to  qualify  its  obvious  merit  by  a  few  key  points  which 
generally  apply  to  all  MM  problems. 

*  The  matrix  dimension  N  (number  of  expansion  modes)  must 
be  taken  sufficiently  large  to  get  a  convergent  result 
but  not  so  large  as  to  cause  numerical  instability. 

*  Small  numerical  errors  due  to  round-off  or  truncation 
usually  have  catastrophic  effects  if  they  cause  apparent 
singularities  or  ill-conditioning  in  the  matrix. 

*  Choice  of  appropriate  basis  functions  not  only  improves 
the  accuracy  but  reduces  the  matrix  size  needed  for  a 
convergent  solution.  In  the  case  of  a  short  dipole,  the 
expected  triangular  current  distribution  suggests  that 
piecewise  linear  or  triangle  basis  functions  should  give 
a  better  representation  of  I(z'). 

*  Galerkin's  technique,  (chosen  expansion  and  weighting 
functions  are  the  same)  is  computationally  convenient  on 
a  computer  but  does  little  to  affect  accuracy.  If 
pulse-pulse  Galerkin  method  is  used  on  the  dipole 
example,  convergence  is  achieved  with  a  different  N. 
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III.  LTSA  MODEL 


In  order  to  examine  the  various  parametric  effects  of 
LTSA  design  on  radiation  patterns,  a  valid  antenna  model  is 
required.  While  LTSA  applications  are  primarily  based  on 
empirical  designs,  recently  published  works  have  advanced  the 
necessary  theory  to  describe  the  LTSA  behavior.  Among  them 
is  a  paper  by  Janaswamy  [Ref.  1].  His  model  enables  the 
prediction  of  pattern  behavior  which  agrees  closely  with  LTSA 
experimental  results.  It  is  Janaswamy' s  MM  model, 
algorithmic  steps  and  Fortran  programs  which  are  used  in  this 
study.  Accordingly,  this  chapter  summarizes  the  purpose  and 
process  of  this  new  model. 

Analogous  to  the  linear  dipole  example,  the  immediate 
objective  is  to  solve  for  the  current  distribution  on  the 
surface  of  the  antenna  for  a  given  excitation.  The  current 
is  determined  by  solving  the  reaction  integral  equation  which 
states  that  [Ref.  1] 

jj  (  J“  •R'r-MB  •H'r)  ds  =  -JJJ  (  Ji*E'r-M1  ‘H^dv  (3.1) 

where  ( E-7- , )  is  the  free  space  field  of  an  electric  test 
source  placed  on  the  surface  S.  (JSM1)  are  the  impressed 
sources  (excitation)  bounded  by  volume  V.  (J^M")  represent 
the  unknown  total  surface  currents  on  the  antenna.  Using  the 
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concept  of  surface  impedance,  the  currents  J“  and  M"  can  be 
related  as 

Ms  =  ZsJs  x  n.  (3.2) 

The  surface  impedance  concept  is  useful  in  simplifying  an 
otherwise  difficult  case  of  the  dielectric  supported  antenna, 
we  recognize  that  Me  =  0  on  perfect  conductors  (air 

dielectric).  However,  on  dielectric  supported  antennas,  the 
presence  of  the  thin  dielectric  slab  is  treated  by  means  of 
the  surface  impedance  Zn.  The  surface  impedance  Z«  is  a 
function  of  the  relative  dielectric  constant  and  the 
substrate  thickness  in  wavelengths  (See  Table  3.1). 
According  to  equation  (3.2),  a  non-zero  value  of  Z«  results 
in  additional  electric  surface  currents.  In  either  case, 
equation  (3.1)  is  solved  for  the  single  unknown 

The  MM  approach  begins  with  selection  of  expansion  modes 
(synonymous  with  basis  functions)  used  to  model  J”.  Besides 
choosing  appropriate  functions,  this  step  has  other 
significant  implications  which  first  demand  consideration. 
The  number  of  modes  should  be  as  small  as  possible  to  reduce 
computation  time  but  large  enough  to  insure  convergent 
results  for  the  radiation  pattern.  To  that  end,  one  helpful 
simplification  is  to  reconsider  the  antenna  structure.  Both 
the  number  of  expanded  unknowns  in 

N 

JB  =  S  InJn,  (3.3) 

n=l 
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TABLE  3.1  SURFACE  IMPEDANCES  FOR  TYPICAL  DIELECTRIC 
SUPPORTED  LTSA  (INDUCTIVE  REACTANCE,  OHMS) 


Substrate 
Thicknes 
d(Ao)  : 

1.04 

Relative 

2.33 

Dielectric 

3.5 

Constant , 

6.0 

er: 

10.5 

.01 

11.86 

11.88 

11.90 

11.94 

12.01 

.  015 

17.82 

17.89 

17.95 

18.09 

18.34 

.021 

25.02 

25.21 

25.39 

25.78 

26.51 

.025 

29.86 

30.19 

30.49 

31 . 16 

32.46 

.030 

35.97 

36.54 

37.08 

38.29 

40.73 

.040 

48.44 

49.84 

51.20 

54.43 

61.65 

.050 

61 . 33 

64 .22 

67.13 

74.55 

94 . 28 
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and  the  computation  time  are  reduced  considerably  when  the 
trapezoidal  plates  of  Figure  3.1a  are  modeled  with  uniformly 
rectangular  patches  as  shown  in  Figure  3.1b.  This  segmen¬ 
tation  allows  overlapping  surface  dipoles  to  represent  each 
Jr,  as  current  flow  along  either  principal  direction  of  the 
plates.  Accordingly,  Figure  3.1c  represents  the  geometric 
approximation  for  each  plate  where  the  feed  gap  is  denoted  as 
2W r.  (This  is  valid  for  flare  angle  2a  less  than  40°. ) 

If  the  number  of  segments  along  H  and  L  are  Nf,  and  Nt., 
respectively,  the  number  of  modes  per  plate  is 

NP  =  Nl(Nh-1)  +  Nh(Nt.-1).  (3.4) 

As  the  two  plates  comprising  the  antenna  surface  are 
identical,  they  both  have  the  same  number  of  modes, 
including  one  mode  to  represent  the  feed  dipole,  the  total 
number  of  modes  on  the  structure  is 

N  =  2Np  +1.  (3.5) 

Selection  of  functions  to  represent  each  mode  aie 
illustrated  in  Figure  3.2  where  the  assumed  direction  of 
current  flow  is  along  the  z  axis.  For  each  surface  dipole, 
the  current  is  represented  by  the  Piecewise  Sinusoid  (PWS)  in 
the  z  direction  and  constant  along  the  y  direction. 

Using  the  representations  of  Figures  3.1  and  3.2,  a  three 
dimensional  view  of  current  is  obtained  by  examining  two 
cross-section  views.  Figure  3.3a  depicts  Cut  AA,  which  is 
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Figure  3.1a  LTSA  Surface  Geometry. 


Figure  3.1b  Surface  Dipoles  Overlapping  on  Rectangular 
Patches. 
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Figure  3.2  Surface  Dipole  i 


Figure  3.3a  Cross  Section  of  Current  Distribution  Along 
the  Length  Dimension  of  the  Plate  (Cut  AA). 


M 


Cut  BB  £=H 


C=o 


Figure  3.3b  Cross  Section  of  Current  Distribution  Along 
the  Height  Dimension  of  the  Plate  (Cut  BB). 
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used  to  portray  current  behavior  along  the  length  dimension. 
Similarly,  Figure  3.3b  is  Cut  BB ,  which  provides  a  view  of 
current  behavior  along  the  height.  Consequently,  current  is 
analyzed  using  two  separate  plots  of  Jtra  and  J  wad  cis  S66ii  in 
Appendix  C. 

The  test  modes  (weighting  functions)  for  J-1-  are  chosen 
identical  to  the  expansion  modes  for  J®,  hence  the  method  is 
Galerkin.  Although  other  choices  are  possible,  PWS  Galerkin 
is  best  suited  in  the  MM  procedure  for  scattering  problems  in 
free  space  and  is  proved  to  be  satisfactory  in  this  model. 

Having  formulated  the  MM  strategy,  equation  (3.1)  is 
solved  by  substituting  (3.2),  (3.3)  and  using  the  familiar 
form 

[Z„„JIJ„]  =  tVm)  (3.6) 


where 


=  Z„  ||  JrWjds  "  |{  J  *  ds  (3.7) 

and 

V"  *  W  -  H^ldv.  (3.8) 

V 

Here,  Zmin  represents  the  mutual  impedance  between  the  mtJl 
testing  mode  (weighting  function)  and  the  nth  expansion  mode 
(basis  function)  and  can  be  termed  the  inner  product, 
<  J„ ,  Ei)i> .  The  complete  discussion  of  equation  (3.7)  need  not 


33 


be  expanded  here  but  two  aspects  of  Zmn  are  significant  in 
this  LTSA  model . 

*  Janaswamy  treats  the  cases  of  parallel  and  non-parallel 
dipoles  with  different  algorithms  and  uses  closed  form 
expressions  where  appropriate.  In  concert  with  the 
problem  symmetry,  his  approach  improves  computational 
efficiency. 

*  Filling  the  impedance  matrix  is  a  time  consuming  task. 
Table  3.2  provides  typical  CPU  times  required  to  fill 
and  perform  LU  decomposition  on  such  a  matrix.  Owing  to 
aforementioned  symmetry,  matrices  require  only  dimension 
N' ,  where  N'  =  Np+l.  Therefore  the  matrix  [Zm„]  has 
( N ' ) 2  elements.  The  preponderance  of  LTSA  modeling 
problems  requires  N'>130,  for  which  a  virtual  memory  of 
1500  Kilobytes  is  recommended.  For  N'  greater  than  200, 
2  Megabytes  of  memory  are  necessary. 

As  suggested  in  the  dipole  example,  increasing  N  improves 
the  approximation  of  the  current  distribution.  However, 
practical  consideration  of  the  required  computer  resources 
dictates  a  limit  for  N.  Janaswamy  has  demonstrated  that 

convergent  results  are  obtained  for  the  radiation  pattern 
when  the  LTSA  height  and  length  are  partitioned  into  four  to 
five  segments  per  wavelength.  This  criteria  was  used 
throughout  the  modeling  study. 

It  is  worth  mentioning  that  an  accurate  LTSA  input 
impedance  can  be  calculated  only  if  convergence  is  achieved 
in  the  solution  for  [In].  Because  an  even  larger  N  is 
required,  this  task  is  not  undertaken. 

The  matrix  [Vm]  in  equation  (3.3)  is  determined  by  the 
type  of  excitation  used  in  the  model.  In  the  present  model, 
a  delta  gap  generator  is  used,  corresponding  to  a 
transmitting  antenna.  By  reciprocity,  the  solution  is  also 
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TABLE  3.2  TYPICAL  CPU  TIME  AND  STORAGE  REQUIREMENTS  FOR  THE 
IMPEDANCE  MATRIX  OF  THE  LTSA  MODEL  (USING  VS 
FORTRAN  AND  IBM  3033) 


Segments 
per  Plate 

Nr.  N„ 

Dimensions 

N 

of  Matrix 

N' 

CPU  Time 
in 

Minutes 

Matrix  Stor¬ 
age  in  Bytes 
( Formatted 
File) 

15 

2 

87 

44 

8.0 

16K 

20 

2 

117 

59 

14.4 

280K 

25 

2 

147 

74 

22.6 

440K 

15 

5 

261 

131 

71.5 

1 . 4M 

23 

6 

495 

248 

245  . 

4 . 9M 
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valid  for  the  receiving  case.  This  choice  results  in  Vm 
equal  to  zero,  except  at  the  source  segment  where  it  is 
unity.  Due  to  the  plate  corners  which  would  cause  undue 
fringing  around  the  feed  gap,  the  normal  direct  connection  is 
not  used.  Instead,  the  generator  is  superimposed  on  a  strip 
dipole  centered  across  the  feed  gap,  as  depicted  in  Figure 
3.4.  The  feed  dipole  length  and  width  are  2Ld  and  2Wd 
respectively.  Results  of  this  study  show  that  the  far-fields 
are  relatively  insensitive  to  changes  in  the  feed  dipole 
dimensions  (See  Figure  3.5).  The  apparent  H-plane  behavior 
change  shown  in  Figure  3.5b  is  attributable  to  the  fact  that 
only  one  expansion  mode  is  used  for  the  source.  Figures  A.l 
through  A. 4  further  demonstrate  more  vividly  that  these 
effects  are  insignificant. 

In  summary  of  the  purpose  and  process  of  the  LTSA  model, 
the  MM  technique  provides  an  approximation  of  the  current 
distribution  which  is  subsequently  used  to  compute  the  far- 
field  radiation  pattern.  A  point  worth  emphasis  is  the 
criterion  of  choosing  four  to  five  segments  per  wavelength  in 
either  principal  direction  of  the  plate.  While  this  is  not 
enough  to  obtain  convergent  results  for  current  distribution, 
it  is  quite  accurate  for  the  radiation  pattern.  Moreover, 
this  model  provides  predicted  patterns  well  corroborated  by 
empirical  results.  Thus,  the  only  apparent  limits  on  its 
effectiveness  are  established  by  the  computer  resources 
available . 
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Figure  3.4  Source  Model  of  Unit  Volt  Excitation 
on  a  Surface  Dipole  Placed  across  the  Feed  Gap. 


37 


H'Plane  8aafflvldth  (Oagraea)  E'Plana  Beamvldth  COagreee) 


Figure  3.5a  E'Plane  Effects  of  Dipole  Modeling. 
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Figure  3.5b  H-Plane  Effects  of  Dipole  Modeling. 
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RoLotlve  Gain,  1st  Sldaloba  (dh)  Rotative  Gain,  1st  Sldelobe  (db) 


IV.  PREDICTED  PATTERN  BEHAVIORS 


A .  BACKGROUND 

The  theoretical  model  provides  the  opportunity  to  develop 
specific  parametric  analysis  of  the  LTSA.  Heretofore,  these 
analyses  have  only  been  experimentally  performed.  Therefore, 
the  ability  to  accurately  predict  radiation  pattern  behavior 
is  a  significant  step  in  the  LTSA  design  methodology. 

Before  the  correlation  of  behavior  changes  to  design 
changes  can  yield  useful  results,  a  validated  set  of  design 
parameters  is  necessary.  For  reference  and  comparison 
purposes,  the  following  LTSA  design  parameters  are  used: 

*  Length,  L  =  3 . OAo 

*  Height,  H  =  0.9Ao 

*  Feed  point  gap,  2W*  =  .02Ao 

*  Flare  angle,  2a  =  12° 

*  Substrate  dielectric  constant,  =  2.33 

*  Substrate  thickness,  d  =  .021Ao. 

The  predicted  radiation  pattern  for  the  reference  data  is 
shown  in  Figure  4.1. 

The  modeled  behavior  corresponds  very  well  to 
experimental  results  obtained  with  this  specific  design  [Ref. 
1J.  Hence,  the  investigation  of  LTSA  behavior  proceeds  by 
varying  one  design  parameter  at  a  time  and  observing  the 
attendant  effects  on  the  radiation  patterns.  This  technique 
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is  well  illustrated  in  Appendix  A  which  contains  radiation 
pattern  comparisons  for  the  designs  modeled  in  this  study. 

A  more  succinct  analysis  is  obtained  by  der  -  th* 

main  beam  and  sidelobe  properties  with  th-  cs 

for  each  plane.  They  are: 

*  Half-power  or  3  db  beamwidths  (HP*.}- 

*  10  db  beamwidth, 

*  Level  of  first  sidelobe  relative  to  main  bean. 

Together  these  attributes  indicate  the  directivity  of  the 
antenna  and  provide  a  basis  with  which  to  compare  antenna 
designs.  The  equation  used  to  approximate  the  directivity  of 
the  LTSA  is 


D  =  52,525/ (HPeHPh)  (4.1) 

where  HPC  and  HPh  are  in  degrees  and  a  smooth  elliptical  main 
beam  cross-section  is  assumed  [Ref.  2:p.  396], 

B.  LENGTH  OF  THE  LTSA 

It  has  been  demonstrated  by  earlier  theory  and 
experiment,  that  for  LTSA  lengths  (L)  of  three  to  ten 
wavelengths,  narrower  beamwidths  are  obtained  by  increasing  L 
(Ref.  4].  Based  on  known  behavior  of  antennas  such  as  the 
Vee  dipole  and  horn  antennas,  this  is  a  reasonable 
expectation.  However,  results  obtained  with  the  present  model 
are  accurate  enough  to  offer  more  specific  observations. 
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Radiation  patterns  for  L  =  3  Ao,  4  Ao  and  5^a  have  been 
computed  for  heights  (H)  of  .  25Ao  to  1.5Ao.  (See  Figures  A. 5 
through  A. 16).  A  synopsis  of  significant  findings  concerning 
length  is  deduced  from  Figure  4.2.  For  a  nominal  height  of 
one  wavelength,  radiation  pattern  characteristics  are  plotted 
as  a  function  of  length. 

As  L  is  increased  the  progressive  decrease  in  beamwidth 
is  evident  in  both  the  E  and  H  planes.  Since  the  directivity 
is  inversely  proportional  to  the  beamwidths'  product,  the 
improvement  between  L  =  3Ao  (HPE  =  35°,  HPM  =  45°)  and  L  = 
5 Ao  (HPe  =  21°,  HPm  =  37°)  is  approximately  3  db . 

Also  noteworthy  is  the  concurrent  sidelobe  reduction 
which  is  most  apparent  in  the  H-plane.  Although  Figure  4.2 
shows  only  the  relative  gain  of  the  first  sidelobe,  review  of 
Figures  A. 11  and  A. 12  indicates  that  the  main  beam  is  well- 
formed  and  that  there  are  no  sidelobe  anomalies.  Suppression 
of  the  E-plane  sidelobes  below  -10  db  is  achieved  with  L  > 

4  Ao . 

Figures  B.l  through  B.6  demonstrate  that  the  observed 
trends  for  changing  L  are  consistent  for  several  heights 
(.25Ao  -  1.5Ao) .  Thus,  a  further  interpretation  regarding 
length  is  worth  consideration.  To  compute  an  actual  antenna 
gain,  the  input  reflection  coefficient  or  the  input  impedance 
is  required.  However,  design  for  a  particular  directivity 
can  be  accomplished  with  equation  (4.1). 
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For  example,  if  a  directivity  of  17  db  is  required,  the 
HPr.  -  HP,,  product  must  be  less  than  1048.  Figure  4.2  indicates 
at  L  =  4^0  that  HP*  =  25°  and  HPW  ~  41°.  It  is  apparent  that 
a  nominal  antenna  length  of  4Ao  is  at  least  a  good  starting 
point  for  this  design. 

C.  LTSA  HEIGHT 

Earlier  theory  for  the  LTSA  had  only  been  developed  for 
large  height  (H  <,  3Ao .  Experiments  with  an  air  dielectric 
LTSA  (L  =  7.2Ao)  showed  that  a  progressively  narrower  E-plane 
beamwidth  is  obtained  when  it  was  shortened  to  1 . 5Ao ,  but 
that  at  . 7 5 Ao ,  the  HP=  broadened  once  again  [Ref.  5].  This 
apparent  resonant  effect  appears  in  modeling  a  shorter 
antenna  (L  =  4  Ao)  as  depicted  in  Figure  4.3a.  When  H  = 
1.25Ao,  there  is  a  distinct  narrowing  of  HPE.  Compared  to  H 
=  1.5  Ao,  this  represents  a  directivity  increase  of  about  2 
db.  For  H  decreasing  below  lAo,  the  HPE  tapers  to  less  than 
20°  at  H  =  . 25Ao . 

In  the  H-plane  the  cited  experiments  showed  a  broadening 
of  HPh  for  decreasing  height.  This  is  consistent  with  the 
predicted  behavior  shown  in  Figure  4.3b.  However,  as  H  is 
reduced  below  .75Ao,  HPh  decreases  and  the  10  db  beamwidth 
virtually  disappears.  Figure  A. 22  shows  the  emergence  of 
large  sidelobes  at  H  =  .25Ao,  while  figure  A. 23  demonstrates 
that  the  beam  is  not  well  formed  at  H  =  . 50Ao-  These  effects 
are  also  observed  at  L  =  3Ao  (Figure  B.7)  and  L  -  5Ao  (Figure 
B  .  9  )  . 
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In  essence,  both  theory  and  experiment  suggest  that  an 
optimum  H/Ao  may  exist  in  a  range  which  depends  on  L/Ao.  For 
example,  Figures  A. 20  and  A. 21  show  that  a  range  of  optimum 
height  appears  very  narrow  -  about  lAo.  On  the  other  hand, 
for  L  =  5  Ao,  Figures  A. 30  and  A. 31  indicate  a  relatively 
broad  range  of  optimum  H  between  1 . OAo  and  1.5Ao. 

It  is  necessary  to  qualify  implications  of  a  so-called 
optimum.  In  this  study,  dimension  multiples  of  A°/4  are  used 
for  L  and  H.  Consequently,  the  interpolation  of  data  points 
must  be  judicious.  It  is  especially  significant  as  the 
relative  contributions  of  other  fixed  design  parameters  may 
change  the  radiation  pattern  with  synergistic  effect. 

Due  to  the  noted  pattern  behavior  attributed  to  changes 
in  height  and  length,  a  possible  relationship  to  the 
length/height  ratio  (L/H)  is  explored  in  Figure  4.4.  No 
discernible  affirmation  or  refutation  is  immediately  evident. 
Consider  L/H  =  4  for  the  three  modeled  antenna  lengths.  The 
characteristics  of  these  designs  are  shown  in  Table  4.1 
These  comparative  results  generally  show  that  a  given  L/H 
ratio  does  not  directly  determine  the  radiation  patterns.  It 
is  reasonable  to  conclude  that  L/  Ao  and  H/Ao  independently 
contribute  to  the  LTSA  behavior. 

D .  FLARE  ANGLE 

The  effects  of  changing  LTSA  flare  angle  (2a)  are 
experienced  in  two  competing  phenomena.  Consider  first  that 
increasing  the  flare  angle  in  similar  traveling  wave  antennas 
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Figure  4.4a  E“Plane  Effects  of  Changing  the  L/H  Ratio. 


6.0  11.0  16.0 
Length  to  Height  Ratio 


Figure  4.4b  H-Plane  Effects  of  Changing  the  L/H  Ratio. 
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TABLE  4.1  PATTERN  BEHAVIOR  FOR  DESIGNS  WITH  L/H  =  4. 
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has  demonstrated  that  HPB  is  proportional  to  1/2 a  and  that  HPh  is 
generally  insensitive  to  flare  angle  changes  [Ref.  6].  These 
results  are  justified  in  part  by  analogy  to  aperture  antennas  in 
which  narrower  beamwidths  are  obtained  by  increasing  the  aperture 
size.  The  second  effect  for  increasing  2a  is  a  broader  beamwidth 
which  results  from  an  increased  phase  velocity.  This  is 
explained  by  noting  that  the  fields  between  the  plates  are 
weakened  as  the  plates  move  farther  apart.  The  phase  velocity 
relates  inversely  to  the  energy  stored  in  the  fields.  Therefore, 
as  the  flare  angle  increases,  the  slow-wave  effect  supporting  a 
narrow  beam  is  reduced  and  the  beam  widens. 

The  net  effects  of  these  two  dynamic  factors  are  summarized 
in  Figure  4.5.  It  is  evident  that  below  21°  the  beamwidths  are 
insensitive  to  flare  angle.  This  contrasts  with  at  least  two 
separate  experiments  for  the  fin  line  LTSA  which  showed  HPh 
unaffected  by  a,  but  also  proved  that  HPE  is  proportional  to  l/2a 
[Refs.  6  and  7].  Another  departure  from  empirical  studies  is  the 
absence  of  beam  symmetry  in  the  predicted  patterns.  It  should  be 
noted  that  in  experiments,  the  beam  symmetry  appeared  in  narrow 
ranges  of  flare  angle  and  only  for  longer  antenna  lengths.  In 
general,  the  relative  contributions  of  the  cited  competing 
factors  seem  to  offset  one  another  in  modeling  below  21°. 

Above  21°,  HP*  increases,  possibly  due  to  the  phase  velocity 
effect  dominating  the  illumination  mechanism.  In  the  H-plane 
there  is  a  concurrent  decrease  in  HPh.  This  may  be  attributed 
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to  larger  current  components  which  move  perpendicular  to  the  H- 
plane  as  the  flare  angle  increases.  Figures  C.l  and  C.23-C.30 
show  the  modeled  current  distributions  for  respective  flare 
angles.  Using  Janaswamy's  model,  surface  impedance  does  not 
change  and  the  current  distribution  appears  to  change  only 
slightly  as  a  function  of  flare  angle.  However,  to  resolve  this 
question  in  further  detail  requires  more  accurate  modeling  of  the 
current  distribution. 

As  a  matter  for  design  consideration,  the  radiation  patterns 
for  different  flare  angles  are  plotted  in  Figures  A.32-A.39. 
Although  these  designs  have  approximately  equal  directivity,  the 
high  sidelobes  are  not  necessarily  acceptable.  For  this  reason 
(as  Figure  4.5  also  indicates)  a  range  of  flare  angles  between 
19°  and  22°  appears  to  be  best,  given  the  same  set  of  other 
design  parameters. 

E .  SUBSTRATE  FACTORS 

As  discussed  previously,  the  plates  are  treated  as  imperfect 
conductors  in  the  dielectric  supported  antenna.  The  surface 
impedance  is  a  function  of  dielectric  thickness  (d)  and  the 
relative  permittivity  (er).  A  modest  range  of  dielectric  designs 
were  modeled  and  the  respective  predicted  radiation  patterns  are 
presented  in  Figures  A.40-A.48. 

It  is  first  noted  that  increasing  "d"  noticeably  reduces  the 
large  E-plane  backlobe  and  slightly  increases  the  relative  gain 
of  the  H-plane  backlobe.  This  is  a  point  of  interest  because 
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backlobes  are  a  minor  disagreement  between  the  LTSA  experimental 
and  theoretical  models.  The  predicted  results  suggest  that  the 
discrepancy  is  due  to  the  manner  in  which  Janaswamy's  model 
approximates  the  effects  of  the  dielectric  substrate  [Ref.  1], 

The  LTSA  behavior  as  a  function  of  substrate  thickness  is 
characterized  in  Figure  4.6.  It  is  observed  that  the  first 
sidelobes  increase  slightly  with  d.  The  E-plane  first  sidelobe 
is  relatively  large. 

The  primary  effect  of  increasing  d  is  the  HPh  improvement. 
Between  d  =  .  OlAo  and  . 05Ao,  HPh  decreases  by  more  than  30%.  The 
E-plane  beam  is  equally  well-formed  but  it  appears  unaffected 
except  where  it  suddenly  narrows  by  30%  between  d  =  .04Ao  and 
.05Ao.  Approximately  2.2  db  improvement  in  directivity  is 
achieved  by  increasing  d  from  ,02lAo  to  ,05Ao. 

The  noted  behavior  is  supported  by  theory  which  allows  that  a 
dielectric  material  typically  slows  the  currents  in  an  endfire 
producing  surface  wave  antenna.  As  the  phase  velocity  decreases 
a  corresponding  narrower  beam  width  is  produced.  The  concurrent 
increase  in  current  is  noted  by  comparisons  among  Figure  C.l  and 
C.31-C.36,  Coincident  with  the  cited  HPE  decrease,  Figures  C.35c 
and  C.36c  show  that  the  transmission  mode  current  at  the  outer 
edge  along  H  does  not  taper  off  as  noted  in  most  other  modeled 
designs.  Because  the  calculated  current  distribution  is  not 
accurate,  it  cannot  be  ascertained  if  this  is  a  meaningful 
observation  or  a  numerical  aberration. 
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The  other  substrate  parameter  to  consider  is  the  relative 
dielectric  constant  (er).  Figure  4.7  depicts  a  proven 

insensitivity  of  pattern  behavior  to  changes  in  er  between  2.33 
and  10.5.  (For  detailed  corroboration,  see  Figures  A.46-A.48  and 
C.37-C.39.)  Simply  stated,  this  occurs  because  the  surface 
impedance  only  depends  on  d  when  the  substrate  is  electrically 
thin,  i.e.,  when  (£r-l)d/  Ao  S  .2.  In  other  words,  for  small  d, 
the  surface  impedance  equation  simplifies  as 

v°  n° 

Z„  =  -  tan{  kodv'Er)  ~  -  k0d  -JTZ  =  70kd  (4.2) 

\/£r  yj £  r* 

This  finding  is  commercially  relevant  because  the  cost  of 
substrate  materials  is  a  factor  in  building  these  relatively 
inexpensive  antennas. 

F.  SUMMARY  OF  RADIATION  PATTERN  EFFECTS 

The  general  purpose  of  modeling  the  LTSA  is  to  gain  insight 
and  apply  the  necessary  theory  to  understand  its  behavior. 
Equally  significant,  however,  is  the  specific  purpose,  which  is 
to  provide  useful  guidance  for  the  design  and  application  of  this 
antenna.  Based  on  the  credibility  of  the  theoretical  model, 
highlights  of  the  parametric  analysis  deserve  reiteration.  Data 
already  exists  for  studies  of  LTSA's  with  large  lengths.  Modeled 
herein  are  design  lengths  of  3Ao,  4Ao  and  5Ao.  The  results  show 
that  directivity  is  directly  and  consistently  proportional  to  the 
antenna  length.  The  design  process  may  commence  by  selecting  an 
appropriate  L  for  a  desired  directive  gain. 
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rlgure  4,7a  E-Plan©  Effects  of  Changing  Substrate  PermlttlvLti 
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Figure  4.7b  H-Plane  Effects  of  Changing  Substrate  Permittivity. 
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The  applied  model  uniquely  provides  an  aid  to  investigate, 
the  height  factor.  A  key  finding  is  that  the  empirically 
observed  optimum  range  of  H/Ao  is  also  evident  in  the  predicted 
pattern  behavior.  Results  further  show  that  while  the  range  of 
an  optimum  H/Ao  does  depend  on  antenna  length,  the  pattern  itself 
depends  on  L/Ao  and  H/Ao  independently. 

For  the  modeled  designs  of  flare  angle  between  6°  and  30°, 
two  mechanisms  appear  to  cancel  the  effects  of  one  another. 
Consequently,  the  directivity  is  virtually  unaffected  by  the 
flare  angle.  Given  other  design  parameters  the  same,  flare 
angles  chosen  between  19°  and  22°  give  the  best  sidelobe 
suppression . 

Substrate  thickness  is  the  dominant  contribution  to  surface 
impedance.  For  electrically  thin  substrate,  er  has  no  effect. 
In  the  modeled  range  of  d,  the  H-plane  beam  narrows  about  5°  per 
.01  Ao.  In  the  otherwise  unaffected  E-plane,  the  backlobe  is 
suppressed  below  -10  db  for  d  l  .  03Ao  and  the  main  beam  suddenly 
narrows  between  d  =  .  04Ao  and  .05  Ao.  Although  beam  sharpening 
can  be  done  with  thicker  substrate,  it  may  be  relatively  costly 
and  impractical. 
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V.  CONCLUSIONS  OF  MODELING  STUDY 


The  scope  of  any  analysis  includes  a  purpose,  supporting 
theory  and  a  demonstrative  application  or  proof.  For  the 
LTSA ,  the  purpose  is  to  develop  a  model  which  accurately 
predicts  its  behavior.  Since  experiments  and  applications 
began  in  1982,  such  a  model  has  become  an  important  goal  in 
the  antenna  community.  The  supporting  theory  advanced  here 
is  to  use  the  MM  to  solve  the  reaction  integral  equation. 
Thereafter,  the  far  fields  are  computed  once  the  current 
distribution  is  known.  The  primary  limitation  to  exploiting 
this  model  is  the  availability  of  CPU  time  required  for 
rigorous  MM  solutions.  Using  this  model,  the  contribution  of 
this  study  is  a  parametric  analysis  of  39  LTSA  designs.  The 
resulting  modeled  behavior  generally  agrees  with  known 
empirical  data.  Consequently,  the  analysis  not  only  is 
consistent  with  applied  theory  but  also  is  useful  for  the 
LTSA  design  process. 

The  presentation  of  ideas  for  review  begins  with  the  MM 
technique.  As  detailed  in  Chapter  2,  MM  makes  it  possible 
and  practical  to  solve  for  an  unknown  function,  f(x),  in  an 
equation  which  does  not  have  a  closed  form  solution.  The 
original  equation  is  essentially  transformed  to  an  equivalent 
matrix  equation  which  is  then  solved  by  conventional  methods. 
There  are  two  major  considerations  in  MM.  Foremost  is  the 
choice  of  basis  functions,  which  ideally  are  selected  to 


faithfully  represent  the  true  f(x).  However,  since  the  form 
of  f(x)  is  often  unknown,  selecting  basis  functions  in 
electromagnetic  applications  is  usually  a  result  of  trial  and 
error.  This  alludes  to  the  second  key  aspect  of  MM.  CPU 
time  may  be  scarce  or  expensive.  the  two  factors  affecting 
the  economy  of  those  resources  are  rate  of  convergence  and 
desired  precision  in  the  solution.  A  good  (or  fortuitous) 
choice  of  basis  functions  improves  the  rate  of  convergence. 
The  matrix  size  (hence  CPU  cost)  grows  geometrically  with  the 
extent  of  required  accuracy. 

The  present  model  employs  the  popular  Piecewise 
Sinusoidal-Galerkin  technique.  The  main  reason  for  this 
choice  is  that  it  helps  to  reduce  the  CPU  time  involved  in 
computing  the  elements  of  the  impedance  matrix.  Furthermore, 
it  has  the  potential  for  reducing  the  number  of  unknowns  for 
a  given  problem.  A  convergent  solution  for  the  radiation 
pattern  (far-field)  is  obtained  with  four  to  five  segments 
per  wavelength.  An  accurate  solution  for  current 
distribution  (near-field)  entails  much  more  CPU  time  because 
more  segments  are  required  for  convergence.  The  model  is 
rigorous  for  the  air  dielectric  case  and  approximate  if  the 
LTSA  has  an  electrically  thin  dielectric  substrate.  To 
precisely  include  the  dielectric  in  the  algorithm  contributes 
little  more  to  the  solution  than  increased  CPU  time. 

In  summation  of  the  parametric  analysis,  there  are 
several  conclusions  which  are  valid  for  the  ranges  of  modeled 
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design  parameters.  Moreover,  it  is  reasonable  to  extrapolate 
behavior  based  on  these  significant  modeling  results: 

*  An  appropriate  antenna  length  is  best  suggested  by  the 
initial  determination  of  required  directivity. 

*  For  a  given  length  there  is  a  range  of  optimum  height 
which  will  contribute  to  a  well  formed  narrow  beam. 

*  Although  the  flare  angle  has  no  effect  on  directivity  in 
a  given  design,  there  is  a  flare  angle  bracket  within 
which  side-lobes  are  ameliorated. 

*  Substrate  thickness  has  a  significant  and  inversely 
proportional  relationship  to  H-plane  beamwidth. 

*  For  electrically  thin  substrate,  the  relative  dielectric 
constant  has  no  effect  on  the  radiation  pattern. 

Future  design  and  development  of  the  LTSA  will  be 
enhanced  by  this  research.  Having  successfully  modeled  the 
pattern  behavior,  the  next  step  is  to  obtain  a  convergent 
solution  for  the  LTSA  current.  Such  a  goal  will  probably 
await  lower  cost  microprocessing  when  improved  accuracy  of  MM 
techniques  is  more  cost  effective.  In  the  interim,  it 

appears  that  the  concept,  the  ideas  and  the  algorithm  used  in 
this  study  may  be  equally  effective  in  the  analysis  of 
similar  antennas. 
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APPENDIX  A  LTSA  RADIATION  PATTERNS 


The  parametric  analysis  of  radiation  patterns  begins  with 
the  known  empirical  design  which  was  accurately  modeled 
(Figure  A.l)  using  the  following  specifications: 

*  Feed  Gap  =  .  02Ao 

*  Length  =  3Ao 

*  Height  =  . 9Ao 

*  Flare  Angle  =  12° 

*  Substrate  Thickness  =  . 02lAo 

*  Substrate  Relative  Dielectric  Constant  =  2.33 

The  study  of  these  parameters  and  their  effects  is  based  on 
Figures  A.2-A.48  which  show  the  computed  radiation  pattern 
for  each  case.  Each  computed  pattern  is  templated  with  the 
reference  pattern  to  illustrate  the  immediate  effects  of  the 
design  change. 

Since  the  theoretical  model  uses  a  strip  dipole  source, 
the  feed  dipole  length  is  varied  as  a  design  parameter  but 
the  feed  gap  is  fixed  at  .02lAo.  The  comparative  effects  of 
this  design  parameter  (Figures  A.2-A.4)  serve  to  help 
understand  the  source  model  rather  than  the  antenna  itself. 

Figures  A.5-A.16  show  calculated  patterns  for  three 
different  antenna  lengths  of  3Ko,  4Ao  and  5Ao.  In  each  case, 
the  height  of  the  antenna  was  changed  in  the  range  .25Ao- 
1 . 5 Ao .  These  same  designs  are  given  a  different  perspective 
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in  Figures  A.17-A.31  which  show  the  effects  of  changing 
height  for  given  antenna  lengths. 

The  range  of  flare  angles  which  were  modeled  is  6°  to 
30°.  Each  modeled  design  is  compared  with  the  original 
design  (12°)  in  Figures  A.32-A.39. 

The  surface  impedance  is  determined  by  both  the  substrate 
thickness  and  permittivity.  Varying  these  parameters  results 
in  radiation  patterns  shown  in  Figures  A.40-A .48. 
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Relative  Power  (db)  Relative  Power 


-180  -ISO  -120  -80  -60  -30  0  33  60  80  120  150  100 

Rngle  Off  Borealght  [Degress) 


Figure  fl.l  Predicted  Patterns  for  a  Known  Empirical  Design. 


-160  -150  -120  -80  -60  -33  0  X  SO  80  120  150  160 

Rngle  Off  Borealght  (Degress) 

Figure  R.2  Effects  of  Decreasing  Feed  Dipole  Length  to  .05X0. 
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-180  -150  -120  -90  -60  -  30  0  30  60  80  120  150  180 

Angle  Off  BoreeLght  [Degrees) 


Figure  A. 3  Effects  of  Increasing  Feed  DLpolB  Lenght  to  .20X0. 


-180  -150  -120  -90  -60  -30  0  30  60  90  120  ISO  180 

Angle  Off  Boreelght  (Degrees) 


Figure  A. 4  Effects  of  Increasing  Feed  Dipole  Length  to  .30Xp. 
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Relative  Power  (db)  Relative  Power  (db) 


“100  -ISO  -120  -90  -60  -30  0  30  60  90  120  ISO  100 

Angle  Off  Boreelght  (Degrees) 


Figure  R.5  Effects  of  Decreasing  Length  to  3^0  for  H-.25XQ. 


-160  -ISO  -120  -80  -60  -30  0  30  60  80  120  150  IBO 

Angle  Off  Boreslght  (Degreee) 

Figure  R.6  Effects  of  Increasing  Length  to  5X„  for  H-.25xo. 
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ISO  “ISO  -120  -SO  -60  -  30  0  30  SO  80  120  150  100 

Angle  Off  Boreelght  (Degree*) 


Figure  FI. 7  Effects  of  Decreasing  Length  to  3X„  for  H-.50Xo. 


Figure  fl.8  Effects  of  Increasing  Length  to  5X„  for  H-.50X0. 
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Relative  Power  (db!  Relative  Power  (db) 


-180  -150  -120  -90  -60  -30  0  30  60  BO  120  150  180 

Rngle  Off  Boreaught  (Degree®) 


Figure  R.9  Effects  of  Decreasing  Length  to  3XD  for  H“.75X0. 


-180  -150  -120  -SO  -to  -X  0  X  $0  90  120  150  180 

Rngle  Off  Boreelght  (Degreee) 


Figure  fl.10  Effects  of  Increasing  Length  to  5X0  for  H-.75X0. 
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Relative  Power  (db)  _  Relative  Power  (dbJ 


-180  -150  -120  -80  -to  -30  0  30  SO  80  120  150  100 

Rngle  Off  Boreslght  [Degrees) 


'Igure  R.ll  Effects  of  Decreasing  Length  to  3X„  for  H“1.0Xo. 


•180  -ISO  -120  -SO  -80  -X  0  X  80  X  IX  ISO  IX 

Angle  Off  Boreslght  (Degrees) 

Figure  FI.  12  Effects  of  Increasing  Length  to  5\0  for  H-1.0Xo. 
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Rotative  Power  (db)  -n  Relative  Power  (dbl 


-180  -150  -120  -80  -60  -30  0  30  8D  80  120  ISO  180 

Angle  Off  Boroslght  (Degrees) 


'Igure  fl.13  Effects  of  Decreasing  Length  to  3XD  for  H“1.25X0. 


-180  -150  -120  -80  -60  -30  0  30  60  80  120  150  180 

Rngle  Off  Boreelght  (Degrees) 


Figure  fl.H  Effects  of  Increasing  Length  to  5X0  for  H-1.25X0. 
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-1B0  -ISO  -120  -80  -60  -30  0  30  60  80  120  ISO  180 

Angle  Off  Boreelght  (Degrees) 


Figure  R.15  Effects  of  Decreasing  Length  to  3X„  for  H-1.50X0. 


-160  -ISO  -120  -80  -60  -30  0  30  60  80  120  150  ISO 

Angle  Off  Boreelght  (Degreee) 


Figure  fl.16  Effects  of  Increasing  Length  to  5X0  for  H-1.50X,. 
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-190  -150  -120  -80  -60  -30  0  30  SO  80  120  ISO  II 

Angle  Off  Boreslght  [Degress) 

Figure  R.l 7  Effects  of  Decreasing  Height  to  .  25Xa  for  L“3X0. 


•180  -ISO  -120  -90  -60  -30  0  30  60  90  120  150  180 

Angle  Off  Sorestght  (Degrees) 


Figure  FI.  19  Effects  of  Decreasing  Height  to  .75Xa  for  L-3X0. 


-100  -ISO  -120  -90  -60  -X  0  X  60  90  120  150  1® 

Angle  Off  Boreelght  (Degreee) 


Figure  FI. 20  Effects  of  Increasing  Height  to  1.25X0  for  L-3X0. 
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Relative  Power  (db)  Relative  Power  (dbl 


-180  “193  -120  -SO  -80  -X  0  X  SO  TO  120  ISO  180 

Angle  Off  Boreelght  [Degrees) 


"Igure  R.23  Effects  of  Decreasing  Height  to  .  50X0  for  L“4X0. 


-180  -isn  -120  -SO  -€0  -X  0  X  60  SO  IX  150  IX 
Angle  Off  Boreslght  (Degrees) 


Figure  R.24  Effects  of  Decreasing  Height  to  .75X„  for  L-4X0. 
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Relative  Power  (db)  Relative  Power  (db) 


-1BO  -150  -120  -80  -60  -30  0  30  60  80  120  ISO  100 

Angle  Off  Boreslght  (Degrees) 


'Igure  R.25  Effects  of  Increasing  Height  to  1.25X0  for  L-4X0- 


-180  -ISO  -120  -SO  -60  -30  0  30  60  BO  120  ISO  100 

Rngle  Off  Boree Ight  (Degrees) 


Figure  R.26  Effects  of  Increasing  Height  to  1.50X0  for  L“4X„. 
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Relative  Power  (db)  _  Relative  Power  (db) 


-180  -iso  -120  -90  -60  -30  0  30  SO  90  120  ISO  100 

Angle  Dff  Borealght  [Degrees) 


'Igure  R.2 7  Effects  of  Decreasing  Height  to  .25Xn  for  L-5X0. 


-100  -ISO  -120  -80  -60  -  30  0  30  60  80  120  150  100 

Rngle  Off  Boreslght  (Degress) 


Figure  R.28  Effects  of  Decreasing  Height  to  .5OX0  for  L“5X„. 


Rotate v®  Power  t db)  _  Rotate vo  Power  (dbl 


"egure  fl.29  Effects  of  Decreasing  Height  to  ,75X0  for  L-5X0. 


Figure  fl.30  Effects  of  Increasing  Height  to  1.25\0  for  L-5xo. 
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Relative  Power  (db)  Relative  Power  tdb) 


“180  -150  -120  -90  -60  -30  0  30  BO  90  120  150  100 

Angle  Off  Boreelght  [Degrees) 


Igure  R.31  Effects  of  Increasing  Height  to  1.50XD  for  L--5X0. 


100  -ISO  -120  -90  -60  -30  0  30  60  90  120  150  100 

Angle  Off  Boreelght  (Degrees) 


Figure  fl.32  Effects  of  Decreasing  Flare  Rngle  to  6°. 
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Relative  Power  (db)  Relative  Power  fdb) 


-180  -ISO  -120  -90  -60  -30  0  30  60  BO  1ZD  150  100 

Angle  Off  Boreelght  [Degrees) 


Figure  fl.33  Effects  of  Decreasing  Flare  Angle  to  9°. 


■180  -ISO  -120  -80  -60  -30  0  30  60  90  120  150  100 

Angle  Off  Boreelght  (Degrees) 


Figure  fl.34  Effects  of  Increi  sing  Flare  flngle  to  15°. 
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Relative  .  over  (db)  Relative  Power  Idb) 


-180  -ISO  -120  -90  -60  -30  0  30  60  BO  120  150  180 

Angle  Off  Boreelght  IDegreee) 


Figure  fl.35  Effects  of  Increasing  Flare  Rngle  to  18°. 
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Relative  Power  (dbl  Relative  Power  Idb) 


-180  -ISO  -120  -90  -60  -30  0  30  60  80  120  150  180 

Angle  Off  Boreelght  (Degrees) 


Figure  fl.37 


Effects  of  Increasing  Flare 


Rngle  to  24°. 


-180  -ISO  -120  -B0  -60  -30  0  30  60  80  120  150  100 

Rngle  Off  Boreelght  (Degrees) 


Figure  fl.38  Effects  of  Increasing  Flare  Rngle  to  27° 
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Relative  Power  (db)  Relative  Power  (dbl 


-180  -ISO  -120  -80  -60  -30  0  30  60  BO  12D  ISO  180 

Rnglo  Off  Boreelght  (Degrees) 


Figure  fl.39 


EFfecls  of  Increasing  Flare 


Rngle  to  30°. 


-180  -150  -120  -BO  -60  -30  0  30  60  BO  120  150  180 

Angle  Off  Boreelght  (Degrees) 


Figure  R.40  Effects  of  Decreasing  Substrate  Thickness  to  . 0 1 0 Xg. 
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Relate va  Power  (db)  __  Relative  Power 


ISO  -ISO  -120  -SO  -60  -30  0  30  SO  BO  120  ISO  ISO 

Angle  Off  Boreslght  (Degrees) 


'Igure  fl.41  Effects  of  Decreasing  Substrate  Thickness  to  ,015X0. 


-ISO  -ISO  -120  -SO  -60  -X  0  X  60  SO  120  ISO  IX 
Angle  Off  Boreelght  (Degrees) 


Figure  FI. 42  Effects  of  Decreasing  Substrate  Thickness  to  .025X0. 
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Relative  Power  (dbl  Relative  Power  (dbl 


“180  -150  -120  -SO  -60  -30  0  30  60  80  120  ISO  100 

Angle  Off  Boreslght  (Degrees) 


'Igure  R.43 


Effects  of  Decreasing  Substrate  Thickness  to  ,D30X0 


-100  -ISO  -120  -80  -60  -X  0  30  60  80  120  150  100 

Angle  Off  Bareelght  (Degrees) 


Figure  fl.44  Effects  of  Decreasing  Substrate  Thickness  to  .040X0 
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Relative  Power  (db)  -n  Relative  Power  tdbl 


Thtckne98  (d)  -  .Q50X„  / 

Thickness  (d)  “  .021^.  / 
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H-Plane 
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‘J  i  i  A  / ' 


E“Plane 
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-160  -ISO  -120  -90  -fio  -X  0  30  ED  SO  120  150  180 

Angle  Off  Boreslght  (Degrees) 


’Igure  FI. 45  Effects  of  Decreasing  Substrate  Thickness  to  .050Xt 


Cn  of  Subetrote  -  3.50 
of_Subetrate  -  2._33 


H-Plane 


-180  -150  -120  -90  -60  -30  0  30  60  90  120  150  1W 

Angle  Off  Boreelghl  (Degrees) 
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Figure  R.46  Effects  of  Increasing  Substrate  Permittivity:  Gp~3.5. 
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Relative  Power  tdb)  Relative  Power  (db) 


-180  -ISO  -120  -90  -60  -30  0  30  60  SO  120  ISO  100 

Angle  Off  Boreslght  [Degrees) 


'Igure  R .47  Effects  of  Increasing  Substrate  Permittivity:  Cp-6.0. 


-180  -ISO  -120  -90  -60  -30  0  30  60  80  120  ISO  ISO 

Angle  Off  Boreslght  (Degrees) 


Figure  R.48  Effects  of  Increasing  Substrate  Permittivity:  Cr-10.5. 
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APPENDIX  B  LTSA  RADIATION  CHARACTERISTICS  AS  FUNCTIONS 
OF  THE  LENGTH  AND  HEIGHT 

The  summary  of  pattern  behavior  is  made  concise,  but 
useful,  by  examining  the  beamwidth  and  sidelobe 
characteristics.  Figures  B.1-B.9  present  this  information 
(in  both  the  E-Plane  and  H-Plane)  versus  the  antenna  length 
and  height.  These  figures  demonstrate  the  results  of 
changing  length  for  different  modeled  heights  (Figures  B.l- 
B.6).  Similarly,  Figures  B.1-B.9  complement  the  analysis  by 
illustrating  the  effects  of  changing  height  for  each  modeled 
length . 
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H“Plans  Beamvldth  (Degrees)  E_Plane  BeamwLdth  [Degrees) 


I  . .  ■  T— — . . . .  "  I 

3.0  4.0  5.0 

Normalized  Length  IX,) 


Figure  B.la  E-Plane  Effects  of  Changing  L;  H  -  0.25X„. 


'i  ■  i  i 

3.0  4.0  5.0 

Normalized  Length  (X,) 


Figure  B.lb  H-Plane  Effects  of  Changing  L;  H  -  O.25\0. 
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Relative  Gain,  1st  Sldelobe  (db)  Relative  Gain,  1st  Sidslobe  (db) 


Beamvldth  (Degrees)  E-Plane  Beamvldth  (Oegroee) 
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Rotative  Gain.  1st  Stdetobe  (db)  Rotative  Pain#  let  Stdetobe  (db) 


H“Plane  Boamvldth  (Degrees)  E”Plane  Beamwldth  (Degrees) 


Figure  B.3a  E_Plane  Effects  of  Changing  L;  H  -  0.75X0. 


i  i  i 

3.0  4.0  5.0 

Normalized  Length  (X,) 


Figure  E.3b  H-Plone  Effects  of  Changing  L;  H  ~  0.75X„. 
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Relative  Gain,  1st  Sldelobe  (dbl  Relative  Gain,  1st  Sldelobe  (db) 


H'PLane  Beamvldth  (Oegreea)  E-Plane  Raamvldth  (Degreee) 


LEGEND 

A  -  3db  Beamwldth 
•  “  lOdb  Beamvldth 
■  -  1st  Sldelobe  Level 


3.0  4.0 

Normalized  Length  (X,) 


Figure  B.5a  E-Plane  Effects  of  Changing  L;  H  -  1.25X0 


LEGEND 

A  -  3db  Beamvldth 
•  ■  lOdb  Beamwldth 
■  -  let  Sldelobe  Level 


3.0  4.0 

Normalized  Length  (X,) 


Figure  B.5b  H~Plane  Effects  of  Changing  L;  H  - 


1.25X,. 
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Relative  Gain,  1st  Sldelobe  (db)  Relative  Gain,  let  Sldelobe  tdb) 


H-Plane  Beamwldth  l Degrees) 


1  - -  I 

LEGEND 

-A 

A  -  3db  Beamwldth 
•  -  lOdb  Beamvldth 
■  -  1st  Sldelobe  Level 

1 

3.0 

- I - 

4.0 

Normalized  Length  (X,) 

1  1 
5.0 

Figure  B.6a  E-Plane  Effects  of  Changing  L;  H  -  1.50XD. 


Figure  B.6b  H-Plane  Effects  of  Changing  L;  H  -  1.50X„. 
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Relative  Gain,  let  Sldelobe  (db)  Relative 


H-Plane  Baamvldth  (Oegraas)  E-Plane  Baamvldth  (Degree*) 


Normalized  Height  (X,) 


Figure  B.7b  H~Plane  Effects  of  Changing  H;  L  ~  3X„. 
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Relative  Gain,  let  Sidelobe  (db)  Relative  Gain,  1st  Sldelobe  (db) 


H'Plana  BaamwLdth  (Degreee)  E~PLone  Beomvldth  (Degrees) 


Figure  B.3a  E~PLane  Effects  of  Changing  H;  L  -  4X0. 
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Figure  B.8b  H-Plane  Effects  of  Changing  H;  L  -  4X„. 
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Relative  Gain,  let  Sldelobe  (db)  Relative  Gain,  let  SI delobe  (db) 


H-Plana  Baamwldth  (Degrees)  E"Plane  Bearawldth  (Degrees) 


Figure  B.9a  E-Plane  Effects  of  Changing  H;  L  -  5X0. 


~r  —  '  i 

0.5  i.o  i.s 

Normalized  Height  (X,) 


Figure  B.9b  H_Plane  Effects  of  Changing  H;  L  “  5X#, 
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Relative  Gain,  let  Sldelobe  (db)  Relative  Gain,  let  Sldelobe  (db) 


APPENDIX  C  LTSA  CURRENT  DISTRIBUTION 


A  key  facet  of  this  unique  planar  antenna  is  the  slow 
wave  structure  which  is  achieved  by  the  use  of  the  dielectric 
substrate.  However,  the  travelling  wave  nature  of  the  plate 
currents  is  responsible  for  obtaining  an  end-fire  main  beam. 
This  effect  is  illustrated  in  the  modeled  current 
distribution  for  each  design  used  in  this  thesis. 

The  current  distribution  plots  in  Figures  C.1-C.39 
provide  a  means  to  correlate  observed  pattern  changes  to  the 
respective  changes  in  current  on  the  antenna.  These 
computations  are  based  on  fair  approximations  sufficient  to 
accurately  predict  the  far-field  pattern. 
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Currant  Ole trlbutLona  on  CUT  BB  (dl  dRflQ  Curront  Cl a tri but Iona  on  CUT  BB 

Igure  C.5  Modeled  Currents  When  Height  -  .25XQ  and  Length  —  3XQ. 
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(c)  JfRfl  Current  Dletri.buti.one  on  CUT  BB  Id)  Jrrq  Current  Di.strlbuti.one  on  CUT  BB 

Figure  C.ll  Modeled  Currents  When  Height  **  0.25Xo  and  Length  -  4\0. 


4.8 


(•Beg)  eaoqj 


(•BOQ)  eaDiy 


<M  O 


n  « 


(OUI)  (f)  ®Pn'nu6oy 


[ouij  (pj  epnvju6Dy 


108 


I 


3.2 


(•Bag)  ssoHd 


(•6eQ)  eao4d 


S 

7  in 


I 


Currant  Di.etri.buti.ona  on  CUT  BB 


SEMI  5  D6  H - V - V - Z-r-yL - V - * - !“LrT-W - H  SEMI 


(■B«q)  esi>4d 


(•Beg)  «soiy 


u)  (D  opnviu6DU 


(dui)  in  spnviuBog 


(•COQ)  BSDLld 


( -Seai  BBDiy 


(Otuj  in  epmuBog 


(dui)  m  epriviuBog 


113 


i 


.  Distributions  on  CUT  88  Id)  Jrrq  Currant  Distributions  on  CUT  BB 

C.17  Modeled  Currents  When  Height  -  0.25Xo  and  Length  -  5X0. 
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CUT  B8  Id)  Current  Distribution®  on  CUT  BB 

Currents  When  HoLght,  —  1.50Xo  and  Length  —  5X0. 
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Lgurs  C.27  Modeled  Currents  When  Flore  Rngle  -  21°. 
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C.30  Modeled  Currents  When 
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Figure  C.31  Modeled  Currents  When  Substrate  Thickness  -  .010X0. 
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Current  Distributions  on  CUT  B8 
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:)  Current  Distributions  on  CUT  B0  (d)  JflflQ  Current  Distributions  on  CUT  BB 

Lgure  C.37  Modeled  Currents  When  Substrate  Relative  Dielectric  Constant  Gp  —  3.5 
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